Cosmology With Axionic-quintessence Coupled with Dark Matter 
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We study the possibility of explaining the late time acceleration with an axion field which is 
coupled with the dark matter sector of the energy budget of the Universe. The axion field arises from 
the Ramond-Ramond sector of the Type-IIB string theory. We study the background evolution of 
the Universe as well as the growth of the matter perturbation in the linear regime. We subsequently 
use the observational data from Sn-Ia, BAO measurements, measurements of the Hubble parameter 
as well as the observational data for the growth of the matter perturbation to constrain our model. 
Our results show that coupled axion models are allowed to have larger deviation from cosmological 
constant by the present observational data. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



I. INTRODUCTION 

Cosmological observations [D-Ql have now established 
the fact that our universe is undergoing a late time ac- 
celerating phase on large scales. This has been confirmed 
by a number of precise cosmological observations. One 
generally believes that an unknown component having a 
large negative pressure is driving this acceleration and ef- 
forts are on to explain this acceleration by suitable mod- 
ification of gravity at large cosmological distances. For 
a detailed review on various model building aspects for 
this late time acceleration, we refer the readers to some 
excellent reviews on this subject @. 

If the acceleration is indeed driven by some unknown 
candidate, then a cosmological constant A with an equa- 
tion of state (e.o.s) w — — 1 is the simplest option. But 
the model is plagued by two very important theoretical is- 
sues: the fine tuning problem and the cosmic coincidence 
problem. Moreover, the observational data itself can al- 
low models which are different from a constant A. These 
are termed as quintessence models [f|, and are primar- 
ily built with scalar fields with sufficiently flat potentials 
(Please see Copeland et al. Q for a detailed review on 
scalar field models). Problem with scalar field models 
are two folds: the energy scale involved with such fields 
is ~ 10~ 3 cV, which is not only much less than the en- 
ergy scale of inflation but also much less than that of 
the supersymmetry (SUSY) breaking. Also, in order to 
achieve the slow-roll for this scalar fields ( which in turn 
ensure the negative e.o.s), the mass of these fields have 
to be of the order of 10~ 33 eV, and this should not be 
corrected to higher values due to SUSY breaking. This is 
extremely difficult to achieve given our understanding of 
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the hierarchy problem in the standard model of particle 
physics. 

The above issues are addressed in the context of string 
theory and a robust model for quintessence was proposed 
by Panda et al. Q (from now on PST model) which is 
based on the idea of axion monodromy in the context of 
Type-II B string theory. The quintessence field here is 
the axion field coming from the zero mode of two form 
and four form fields of Ramond-Ramond sector. When 
the shift symmetry associated with the axion fields is ap- 
proximately broken, in a non-trivial way, the field devel- 
ops a linear potential. This is perhaps the first controlled 
model for quintessence which is carefully constructed in 
the context of string theory. Subsequently the cosmolog- 
ical constraint on this model was derived by Gupta et 
al. @, where it was showed that the observational data 
(coming mainly from the background Universe) allow the 
model to behave very close to the ACDM model when 
the initial value of the axion field is appropriately fixed. 
Changing this initial value can make deviation from the 
ACDM model. 

Another challenge for any quintessence model is to keep 
the coupling of these scalar fields to the standard model 
fields sufficiently weak. This is to prevent the long range 
effect that these scalar fields can have ( due to their small 
mass) in our Solar system. In the PST model, this has 
been successfully achieved since the axion field from the 
Ramond-Ramond sector does not interact with the fields 
of the standard model sector. But such an axion field, 
in principle, can interact with dark matter sector of the 
Universe. At present, our understanding of the dark mat- 
ter sector is very limited, specifically its interaction with 
other fields. Currently, the dark matter is assumed to be 
weakly interacting matter particles (WIMP). Thus it is 
plausible to assume that the axion field, in the present 
context, can interact with the dark matter sector. This 
interaction could be also of non-perturbative nature. 

In this paper, we investigate the cosmology with the 
PST-model of quintessence where the axion field is cou- 
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pled with the dark matter sector of the theory but does 
not couple with visible matter sector of the Universe. We 
first study the background evolution in this set up and 
discuss the possible modifications that these coupling can 
result in the evolution of the Universe. We also study 
the evolution of the matter perturbation in linear theory. 
Finally we use different observational data to study the 
constraint on different parameters in our model. 

In the next section, we formulate the equations for the 
background evolution of the Universe assuming the axion 
field is coupled to the dark matter sector of the Universe. 
In section III, we formulate the matter perturbation the 
linear regime in our model. We discuss the evolution of 
the linear growth of the matter fluctuations. In section 
IV, we discuss the observational constraints in our model 
using various observational data. Finally we conclude in 
section V. 



II. BACKGROUND EVOLUTION WITH 
COUPLED AXION 



We start with an interacting picture where the axion 
field is coupled with the dark matter sector (d) of the 
Universe. The visible matter sector (b) is not coupled 
with the axion field. The relevant equations are given by 
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FIG. 1: Behaviour of the equation of state for the axion field. 
Bottom to Top: W = 0,0.01,0.03,0.06. Q d o = 0.23, Q b o = 
0.05, Xi = 0.7. 
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together with the flatness condition 
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Here C(4>) represents coupling parameter between the ax- 
ion field and dark matter. Due to our ignorance about 
the detail physics for this underlying interaction between 
the axion field and the dark matter, we address this is- 
sue at a phenomenological level and assume C(4>) to be 
a constant following the earlier work by Amendola [Ioj |. 
For C = we recover the uncoupled case. 
Next, we define the following dimensionless variables: 
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FIG. 2: Behaviour of the equation of state for the axion field. 
Bottom to Top: A, = 0.05, 0.3, 0.5, 0.6, 0.7. fl d0 = 0.23, Sl b0 = 
0.05, W = 0.06. 



The density parameter 17^ and the equation of state 
for the axion field are given as, 



^0 =x 2 + y 2 
7 =1 + u\£ = 



2x 2 



x z + y z 



(5) 
(6) 



With this, one can form an autonomous system of equa- 
tions: 



% =W^/3^(1 - - z 2 ) + 3^(1 - - 7) 



V s ^ 
-37(2-7) 



z = - -zfl^l - 7) 



A' =v /3^7A 2 (l-r), 



(7) 



where W = ~ ■ We evolve the above system of equations 
from the decoupling era (a = 1CU 3 ) to the present day 
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FIG. 3: Behaviour of the density parameters for different com- 
ponents. fi d0 = 0.23, n b0 = 0.05, Xi = 0.7, W = 0.06. 

(a = 1). Given the initial conditions for 7, f2^ z and 
A, we can solve the system. The scalar field is initially 
frozen due to large Hubble damping, and this fixes the 
initial condition for 7, 7.; rj 0. The initial value for A, 
Xi, is a parameter in our model. Xi determines the devi- 
ation of w,/, from the initial w<p ~ — 1 frozen state as the 
universe evolves with time. For smaller Xi, the deviation 
extremely small, and the scalar field behaves as a cosmo- 
logical constant at all time. For larger values of A^, Wrf, 
has large deviation from —1 as the universe evolves. In 
general, the contribution of scalar field to the total energy 
density is negligible in the early universe. Nevertheless 
one has to fine tune the initial value of f2^ in order to 
have its correct contribution at present. In this regard, 
Qj>(initial) is related to the Q,j,(z = 0). Similarly the ini- 
tial value of z ( which is related to the density parameter 
for visible matter) can be related to the present day value 
of the fib which is fixed at Qbo = 0.05 in our subsequent 
calculations. 

In our study the scalar field is the axion field con- 
structed using axion monodromy in the Type-IIB string 
theory, the form of the potential for this field is given by 

V(cf>) = ^<j>, (8) 

J a 

where 

n 2 M 2 

fa = -Qff- . M = Ml +M2- (9) 
Here Ml = {2 ^g"v aild ^ 2 = CM SB^ A ° (fij)- 9s is 

the string coupling constant and V\ — L 6 a 3 is the volume 
of the internal space in the ten dimensional space time. 
A is related to the warp factor at the location of NS5 
brane and R is the radius of the Ads like throat ( Sec [|| 
for the detail construction of this model). It can be shown 
that the [12 term have the dominant contribution in the 
potential V(<f>). With this form of the potential, r = in 
the system of equations (7). In figure 1 and 2, we show the 



FIG. 4: Behaviour of the growth factor f — d j}? 3St ) as a 
function of the scale factor. From top to bottom, ACDM, 
W = 0, 0.03, 0.06. n d0 = 0.23 and fi M = 0.05. Xi = 0.7 for 
models different from ACDM. 

behaviour of the equation of state for the axion field 
as a function of redshift keeping A^ fixed and varying W 
( figure 1) and keeping W fixed at varying A, ( figure 2). 
From these figures it is clear that deviates more from 
w = —1 for larger Xi as well as for larger W. For smaller 
values of these two parameters, behaves very close to 
the A value w = — 1. In figure 3, we show the behaviour 
of the density parameters for different components as a 
function of scale factor. It is evident that in the early 
era, Universe is mostly dominated by the visible as well 
as dark matter components. But as the Universe evolves, 
the axion field starts dominating in the late epoch. 



We mention that the axion field is initially frozen due to 
large Hubble damping and that sets the initial w^i = — 1. 
But there is sudden jump initially in the from this 
frozen value as can be seen from the figures 1 and 2. This 
is due to the coupling between the axion field and dark 
matter component. But it does not create any pathology 
in Universe's evolution; because during the time when 
this sudden jump happens the contribution from the ax- 
ion field to the total energy budget of the Universe is neg- 
ligible. So the axion field does not affect the Universe's 
evolution during this time. We should point that similar 
things happens when a thawing scalar field is coupled to 
the total matter component (visible + dark matter) of 
the Universe [i~2| . 



III. LINEAR PERTURBATION 

In this section, we write down the equations that governs 
the growth of fluctuations for dark matter and visible 
matter. In our model the dark matter is coupled to the 
axion field whereas the the visible matter is uncoupled. 
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We work in the longitudinal gauge, 

ds 2 = a 2 [-(1 + 2$)dr 2 + (1 - 2$)dx i dx i ] , (10) 

where r is the conformal time and $ and $ are the 
two gravitational potential. It is well known that in the 
absence of any anisotropic stress $ = \&. We follow the 
prescription given by Amendola and write the equa- 
tions for the density perturbations for visible matter and 
dark matter in the Newtonian limit which is valid for 
small scales. In this scale, one can safely assume that the 
axion field docs not cluster otherwise it will behave like a 
massive dark matter. With these assumptions, the equa- 
tions governing the growth of the linearised fluctuations 
in dark matter and visible matter are given by: 

( H! \ 3 

5 d + ( 1 + ~ 2 PdA s'd - ^(idd5 d n d + -f db 5 b n b ) = o 

(11) 

s' b ' + (i + ^ - 2/3 b x^j s' b - ^( ldb s d n d + Jbb s b n b ) = o 

(12) 



Here = 1 + 2j3ij3j and % is conformal Hubble param- 
eter % = aH. For our case, (3 b = and (i d = W. The 
density fluctuations S d and S b are defined as Si — 
The total density fluctuation for the matter component ( 
dark + visible) is given by 



<> d Q d + S b fl b 
Q d + fl b 



(13) 



In figure 4, we show the behaviour of the growth function 
/ = d d li°oga] as a function of scale factor. One can see that 
the growth is suppressed as one increases the coupling 
constant W . 



IV. OBSERVATIONAL CONSTRAINTS 



In this section, we use the current observational data to 
put constraints on our model (see fill ] for earlier works 
on observational constraints for coupled quintessence). 
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FIG. 6: Allowed behaviour for the equation of state of the axion field at 68% and 95% confidence level. Top Left: uncoupled 
case W = 0, Top Right: coupled case W = 0.01, Bottom Left: coupled case W = 0.03, Bottom Right: coupled case W = 0.06. 
Light Shade for 68% and dark shade for 95%. The dashed line represents the best fitted behaviour and the solid line represents 
line represents the phantom divide. 



We start with the Supernova Type-la observations which 
directly probes the cosmological expansion. It actually 
measures the luminosity distance of different supernova 
explosions at different rcdshifts which is defined as: 



d L (z) = (! + «) 



1 



H(z') 



dz'. 



(14) 



The actual observable, distance modulus is given by 

d L 



fi = 5log— h 25. 

Mpc 



(15) 



We consider the Union2.1 compilation containing 580 
data points for [i at different redshifts (l3j . 

Next, we use the observational constraints on Hubble 
parameter which have been recently compiled by Moresco 
et al. [14| in the redshift range < z < 1.75 using the 
differential evolution of the cosmic chronometers. They 
have built a sample of 18 observational data point for 
H (z) spanning almost 10 Gyr of cosmic evolution. 

Next, we use the combined BAO/CMB constraints as 
derived recently by Giostri et al. [ijj]. This gives the 



constraint on the angular scales for the Baryon Acous- 
tic Oscillation peak in the matter power spectrum. At 
present we have measurements from SDSS survey flfjj , 
6dF Galaxy Survey [lU and more recently, by the Wig- 
gleZ team [i~7l |. The corresponding covariance matrix is 
given by Giostri et al (l8j . 

We further use the measurement for the growth pa- 
rameter f = d j°9 s t The list of all the currently available 

J a log a J 

growth data is given in the reference [19| ( see also refer- 
ences therein). 

With this, we calculate the combined likelihood: 

- 2 log C = -2 log(£ s „ x C bao x C hub x C gr ). (16) 



This likelihood function is a function of the model param- 
eters Qdo,\i and W. As we mentioned earlier, we have 
fixed fib = 0.05 for the visible matter. We maximize the 
likelihood function for W — for the uncoupled case and 
W = 0.01,0.03,0.06 for the coupled case and draw the 
corresponding confidence contours in the fldo — Xi plane 
where tt m0 = Q d o + ^bo- 
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In figure 5, we show the confidence contours in the 
Qdo — K plane for different values of the coupling con- 
stant W. We fix fifcQ = 0.05 for our calculations. The 
first thing to be noticed is that, for the uncoupled case 
(W = 0), ttdo and A^ are completely uncorrected. As 
one increases the value of W, the two parameters become 
slightly anti correlated. One concludes from these plots 
that in the coupled case the upper bound on the initial 
slope of the axion potential Ai is comparatively smaller 
than the uncoupled case. As the value of the axion field 
is inversely proportional to the slope of the potential in 
our model, this translates to a bigger lower bound for the 
axion field in the coupled case. 

We have discussed in the earlier section that the devi- 
ation of the cosmological evolution from the concordance 
ACDM behaviour is controlled by the two parameters Ai 
and W. Hence we use the standard error propagation 
technique [2(| to calculate the error in the equation of 
state w,/, knowing the errors in Ai and ftdo for different 
values of W. The results are shown in the figure 6. In 
these figures, we also show the region in the phantom 
side as obtained by the error propagation. But for our 
axion field the relevant region is above the line = — 1 . 
We observe from these plots that for the uncoupled case, 
the allowed behaviour for w^, is very close to the A line 
w = — 1. This is consistent with the earlier results ob- 
tained by Gupta et al. @ . As one increases the coupling 
parameter W, larger deviation from w — — 1 is allowed. 
This confirms that the coupling with dark matter can 
result the axion-quintcsscncc field to behave differently 
from the cosmological constant A. 

V. CONCLUSION 

The PST model proposed Panda et al. is perhaps the 
first controlled model for quintessence, built in the con- 
text of string theory using the idea of axion monodromy. 
In this set up, the axion field arises in the Ramond- 



Ramond sector of string theory avoiding the interaction 
with standard model fields. But such an axion field can, 
in principal, interact with the dark matter sector of the 
Universe. In this work, we study the cosmology in a PST 
model where the axion field is coupled with dark matter 
sector but not with the visible matter sector. Due to the 
lack of understanding about the detail physics for this 
interaction, we assume the coupling function to be a con- 
stant for simplicity. This kind of coupling function has 
been earlier used by Amcndola and others for studying 
the coupled quintessence models. 

In our model, the axion field is initially frozen due to 
large Hubble damping and the equation of state w$ ~ — 1 
initially. But as the Universe evolves, w^, deviates from 
this frozen value. This deviation is larger for larger cou- 
pling and larger initial slope of the potentials A, . Next we 
set up the equation of for linear matter density perturba- 
tion in our model and study the evolution of the growth 
function showing that stronger coupling suppresses the 
growth. 

Finally we use the latest data from Supernova Type-la, 
Baryon acoustic oscillations, measurements of the Hubble 
parameter as well as the measurement of the growth, to 
constrain the model parameters. For the uncoupled case, 
we confirm the earlier results by Gupta et al.Q that data 
allow extremely small deviation from the cosmological 
constant A. But as one introduces the coupling between 
the axion field and the dark matter, larger deviation from 
the w = — 1 behaviour is allowed with increasing strength 
of the coupling. 
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